The longest excursion of stochastic processes in nonequilibrium systems 
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We consider the excursions, i.e. the intervals between consecutive zeros, of stochastic processes 
that arise in a variety of nonequilibrium systems and study the temporal growth of the longest one 
^max(i) up to time t. For smooth processes, we find a universal linear growth (Z ma x(t)) — Qoo t with a 
model dependent amplitude Qoo- In contrast, for non-smooth processes with a persistence exponent 
9, we show that (Z m ax(i)) has a linear growth if 9 < 6 C while (Z m ax(t)) ~ t}~^ if 9 > 9 C . The amplitude 
Qoo and the exponent ip are novel quantities associated to nonequilibrium dynamics. These behaviors 
are obtained by exact analytical calculations for renewal and multiplicative processes and numerical 
simulations for other systems such as the coarsening dynamics in Ising model as well as the diffusion 
equation with random initial conditions. 
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Introduction. Nonequilibrium dynamics in many-body 
systems keep offering new challenges despite several 
decades of research. An example of such a system, 
among others, is the Ising model undergoing phase or- 
dering after a rapid quench in temperature In such 
systems, the relevant stochastic process X(t) that rep- 
resents, at a fixed point in space, the evolving spin in 
the Ising model (or e.g., the field in the diffusion equa- 
tion) is generically a complex one with nontrivial history 
dependence. Traditional two-time correlation function 
(X(ti)X(t,2)) is typically not sufficient to characterize 
the complex temporal history of such a process. One 
simple measure of this history dependence that has at- 
tracted much attention in the recent past, both theo- 
retically 0, Q and experimentally [H , is the persistence 
p(ti,t 2 ) defined as the probability that the process X(t), 
adjusted to have zero mean, has not changed sign in the 
interval \t\,t2\. In several such nonequilibrium systems 
persistence p(to,t), for t ^> to, decays as a power law, 
Po(t) = pftoA 3> to) ~ t~ e , with a nontrivial persistence 
exponent 9 [3(. 

A stochastic process X(t) (depicted schematically in 
Fig. 1) evidently does not change sign between two con- 
secutive zero crossings. The persistence po(t) is simply 
related to the probability distribution of time intervals 
(or excursions) between successive zeros and is clearly 
one, but not the only one, possible characterization of 
the history dependence of X(t). In this Letter we pro- 
pose an alternative yet simply measurable characteristic 
of the history of X(t) via an extreme observable that elu- 
cidates, in a natural way, the important role played by 
extreme value statistics in such nonequilibrium systems. 
In particular, our results illustrate the universal features 
of extreme statistics in generic many-body nonequilib- 
rium systems and provides, in addition, interesting con- 




FIG. 1: Intervals between zero-crossings (excursions) of a 
stochastic process X(t). 



nections with the theory of records that attracted much 
interest recently in the context of random walks 1 5] , grow- 
ing networks [6[ and pinned elastic manifolds [7|. 

For a typical realization of a generic process X(t) with 
N = N(t) zeros in the fixed time interval [0,t] (see Fig. [TJ, 
let {ti,T2, • • • ,tn} denote the interval lengths between 
successive zeros and A(t) denote the length (or age) of 
the last unfinished excursion. Our proposed extreme ob- 
servable is the length of the longest excursion up to t 



imax(i) = max(ri, t 2 , • • • , t n , A(t)) 



(1) 



Let us first summarize our main results. We find 
rather rich universal late time behavior of the average 
(lmnx(t)) for generic stochastic processes X(t) arising in 
nonequilibrium systems. Such processes typically belong 
to two broad classes Q: smooth (i.e., with a finite den- 
sity of zeros) and non-smooth (with infinite density of 
zeros). While persistence typically decays algebraically, 
Po(t) ~ t~ e , irrespective of the smoothness of the pro- 
cess, (Jmax(i))) in contrast, turns out to be sensitive to 
the smoothness of X(t). For smooth processes, (l max (t)} 



2 



always grows linearly with time 

(^max(^)} — Qoo t ; 



(2) 



albeit with a model dependent prefactor > 0. In con- 
trast, for non-smooth processes, it grows either as in ([2]), 
or as 



(imax(t)> ~ t 



(3) 



where the exponent < ip < 1 (sublinear growth), de- 
pending on whether the associated persistence exponent 
9 of the process is less (9 < 9 C ) or greater (9 > 9 C ) than a 
critical value 9 C . We establish these behaviors via exact 
analytical calculations for two simple models, one corre- 
sponding to each class: multiplicative (smooth) and re- 
newal (non-smooth) processes. In addition, we perform 
extensive numerical simulations in a variety of nonequi- 
librium systems, including the diffusion equation with 
random initial conditions and coarsening dynamics of the 
Ising model both below and exactly at the critical tem- 
perature T = T c . For the latter case, our results (0 E|) 
provide new universal quantities associated to nonequi- 
librium critical dynamics. 

The full knowledge of the distribution of i max 8] allows 
in particular the computation of its average. However, 
here we compute (l max (t)) by using the relationship 



d{l max (t))/dt = Q(t) 



(4) 



where Q(t) denotes the probability that the last excursion 
in [0,t], A(t) in Fig. [T] is the longest one 



Q(t) = Prob[; max (i) = A(t)} 



(5) 



Thus Q(i) is the rate at which the "record" length of an 
excursion is broken at time t. Indeed, if the total interval 
increases from t to t + dt, the random variable l max (t) 
either increases by dt (if the last excursion happens to be 
the longest one and the probability for this event is Q(t)) 
or stays the same (with probability 1 — Q(t)). Taking 
average gives ((4|). Henceforth we focus on Q(t), rather 
than (imax(t)) directly. 

Renewal processes. Let X(t) be a renewal process with 
successive interval lengths Tj 's being independent random 
variables, each distributed according to a Levy law with 
parameter 9, p(r) ~ t^ 1 ^ 9 , for large r The per- 
sistence is simply po(t) — J°° drp(r) ~ t~ e . The joint 
distribution 5at(ti, T2, • • ■ , rjy, A[t); t) of the intervals de- 
picted in Fig. [1] is then 

qwin, ■ ■ ■ ,T N ,A(t):t) = p(n)p(T 2 ) ■ ■ ■ p(T N )p Q (A(t)) 
x 5(n + t 2 H h t n + A{t) - t) , (6) 

where the S function ensures that the total interval length 
is t. For the last interval to be the longest, the others 
must be shorter than it implying that Q(t) in Eq. ([5]) is 



Q(t) = 



N=0 



db I dr\. 



dT N q N {ri, ■ 



,r N ,b;t) . 
(7) 



Taking Laplace transform of Eq. ||7J) gives a simple form 
for Q0) = f °° dte~ st Q(t) : 



Q(s) 



db- 



p (b)e 



-sb 



o 1 - J Q drp(T)e- 



(8) 



Using p(r) — ~p'q(t) in an integration by part, followed 
by change of variables b = x/s and r — y/s, lead to an 
expression convenient for late time asymptotic analysis 



Q(s) = - 

S JO 



dx 



p (x/s)e x 



Po(x/s)e x + J dyp Q (y/s)e v 



(9) 



For 9 < 1, one can take the limit s — > directly 
in Eq. ([9]), using po(t) ~ {to/t) e for large t with some 
non-universal microscopic time scale to. Interestingly, to 
cancels between the numerator and the denominator in 
Eq. ©, yielding Q(s) - Q*/s and thus Q(t) -> 
for large t (the superscript R refers to renewal process), 
as announced in Eq. ([2]) with a universal constant that 
depends only on 9 (and not on other details) 



dx 
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(10) 



A special case of this general result, 9 = 1/2, corre- 
sponds to Brownian motion if one considers only "large" 
excursions, i.e. r^'s in Fig. [1] larger than some cut-off 
t c . This recovers in a simple way the result Q^(l/2) = 
0.626508..., derived previously by mathematicians [10( 
using rather complicated, albeit rigorous, method. Note 
that Q^(9) in Eq. (JTDJ> vanishes as 6 -> 1. A plot of 
Q*(9) vs. 9 is shown in Fig. Hi. 

In contrast, for 9 > 1, the naive substitution of po(t = 
y/s) ~ {to/t) e in the integral in the denominator of Eq. 
© is problematic since the integral diverges. Instead 
a careful analysis shows that J Q X dypo{y/s)e~ v oc (r) s 
as s — * where (r) = f °° dTTp(r). This yields, after 
simple algebra, Q(s) - s" 1 / 6 and thus Q(t) - t~ 1+1/e , 
as announced in Eq. ([3]), with 



ip = 1-1/9 



(11) 



Thus for renewal processes the change of behavior of 
(lmax(t)) happens at 9 = 9 C = 1. Qualitatively this tran- 
sition can be understood by simple scaling arguments 
combining extreme value statistics with the behavior of 
the sum of independent Levy variables 0]. The fact 
that the asymptotics of (Z m ax(i)) for Brownian motion, a 
highly "non-smooth" process with infinite density of zero- 
crossings, corresponds to a special case (9 = 1/2) of the 
renewal process suggests that the latter might qualita- 
tively lead to a good approximation of (l ma , x (t)) for other 
non-smooth processes such as the coarsening dynamics 
of the Ising model, and leads us to hypothesise that this 
change of asymptotic behavior of (lm ax (t)) at a certain 9 C 
might be generic for non-smooth processes. Such an ap- 
proximation of the phase ordering of the Ising model by 
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a renewal process is also useful for other observables 11 1 . 
However, for generic non-smooth processes with 9 > 9 C , 
the scaling relation ip = 1 — 1/6 obtained under renewal 
approximation is in general not valid and tp seems to be 
a new exponent. Numerical results indeed support this 
hypothesis. 

Multiplicative processes. A process X(t) is multiplica- 
tive if the locations of its zeros {ti, t2, • • ■ } are such that 
the successive ratios Uk = ifc-i/tfc are independent ran- 
dom variables, each distributed over U <G [0, 1] with den- 
sity p(U), While the calculation of Q(t) is difficult for 
arbitrary p(U), it turns out that for the special family of 
density parametrized by 6, p(U) = 9U e ~ 1 , one can use 
recent results of Ref. to show that Q(t) — > (where 
M refers to multiplicative process), leading to a linear 
growth of (l max (t)) as in Eq. © with 
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ds t 



-6E(s) 



(12) 



where E(s) — J°° dxe x /x. In particular, for uni- 
form distribution, 8 Q^(l) = 0.624329..., the Golomb- 
Dickman constant that also describes the asymptotic lin- 
ear grow th of the longest cycle of a random permuta- 
tion [lj]. In Fig.[2ji we show a plot of Q^(9). At variance 
with renewal processes, (Z max (t)) — Q^(9) t for all 9. 

To appreciate this result in a more general context, we 
note that a multiplicative process X{t) is non-stationary 
by construction. However, when plotted as a function 
of T = ln(t), the process becomes a stationary renewal 
process in T since the successive intervals on the T axis 
Tk — Tk-\ become statistically independent. Similarly, it 
turns out that for many noncquilibrium processes (e.g, 
diffusion equation with random initial condition), the 
original non-stationary process in real time t becomes 
stationary in T = ln(i) 3] and then the renewal ap- 
proximation is precisely equivalent to the Independent 



3j, known to be a very 
I . For such smooth pro- 



Interval Approximation (IIA) 
good one for smooth processes 
cesses then, multiplicative process is a good approxima- 
tion in real time t. Within the IIA, the interval distri- 
bution in log-time T decays as ~ exp[— 8T] for large T 
where 9 is the associated persistence exponent. In real 
time t, this then corresponds to a multiplicative process 
with p(U) ~ 9U e ~ 1 for small U. If one assumes fur- 
ther that this power law form of p(U) holds over the full 
range of U S [0, 1] one arrives precisely at the model 
studied above with the parameter 9 being the persis- 
tence exponent. Thus, the multiplicative process with 
p(U) — 9U e ~ 1 seems to be qualitatively a good repre- 
sentative of generic smooth processes, leading to the hy- 
pothesis of the asymptotic linear growth of {l ma x(t)} for 
such smooth processes. This hypothesis is supported by 
numerical simulations. 

Numerical results. We have computed Q(t) for various 
processes for which the persistence exponent 9 is known 
either exactly or numerically. Guided by our analytical 



results, we have considered both non-smooth and smooth 
processes and the numerical results are consistent with 
the two broad behaviors announced in Eqs. j2|) and ((3]). 
In the first case, Q(t) — > Qoo as in Eq. |[2}, as shown in 
Fig. [2h while in the second case, Q(t) ~ t~v as in Eq. ([3]), 
as shown in Fig. ^Bp. 

As a prototype of non-smooth processes, we have stud- 
ied the magnetization in the coarsening dynamics of a 
d-dimensional ferromagnetic Ising system of linear size L 
consisting of L d spins o~i — ±1, with periodic boundary 
conditions (pbc). Starting from a random initial condi- 
tion, the spins evolve via Glauber dynamics with nearest 
neighbour Ising Hamiltonian -ffising = — Ylu j) (J i (J j- Our 
results are summarized below: 

• In d = 1 and at zero temperature, we have computed 
Q(t) for the local magnetization X(t) = <Ji(t), for which 
9 = 3/8 Q. Fig. Hi shows a plot of Q(t) vs. t for 
L = 128 and 256. These data show that Q(t) — > Qoc 
with Qoo = 0.725(5), which is very close to the analytical 
value obtained for a renewal process in Eq. (flT)]) with 
9 = 3/8, for which Q£(3/8) = 0.726531.. (while for a 
multiplicative process one has Q^(3/8) = 0.80338...), 
see Fig. Efc. 

• We obtained a similar behavior, i.e. Q(t) — ► Qoo 
for the global magnetization M(t) = L~ d J2i< J i(t) both 
in d = 1 at T = (for which 9 = 1/4 0) and in d = 2 
at the critical point T = T C (for which 9 = 0.237(3) [Hj]). 
These results are shown in Fig. Note that for the 
global magnetization, the agreement between the numer- 
ical value of Qoo and the corresponding Q 00 {9) is only 
qualitative. 

• The exact solution obtained for renewal processes 
shows that if 9 is large enough (in that case 9 > 1), Q(t) 
decays to zero as a power law Q{t) ~ t~^ (JH [Tl]). Re- 
cently, it was shown that for critical dynamics of Ising 
systems starting from a completely ordered state, the 
persistence exponent associated to X(t) = M(t) — (M(t)) 
where M(t) is the global magnetization can be large, for 
instance 9 — 1.7(1) in 2d [16(. Unfortunately, the nu- 
merical computation of Q{t) is quite difficult in that case 
because the exponent ip is seemingly positive but very 
small. Alternatively, starting from a fully magnetized 
state, one can instead consider, as in Ref. [17| , the process 
X(t) = M t (t) - (M(t)) where Mj(t) is the magnetization 
of a line, for which the persistence exponent is even larger 
6 ~ 3.3 In Fig. [2p, we show a plot of Q(t) for this 
process for two different system sizes L — 64, 128. This 
plot is compatible with a power law decay Q(t) ~ 
with t/j — 0.34(1), which is actually rather far from the 
value obtained for a renewal process in Eq. (JTrjl which 
gives rp = 1 - 1/0 ~ 0.7. 

These results for non-smooth processes in coarsening 
dynamics, summarized in Fig.[2fc, are qualitatively and in 
some cases even quantitatively (see Fig. [2k) in agreement 
with the results for renewal processes in Eqs. (fTUl [TT]) . 

As a prototype of smooth processes, we have stud- 
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L=128 
L=256 
Q R „(3/8) 
Q M „(3/8) 
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FIG. 2: a): Q(t) as a function of t for the local magnetization 
in the Ising chain evolving with Glauber dynamics at T = 0. 
b): Q(t) as a function of t on a log- log scale for the magne- 
tization of a line in the 2d-Ising model at T c evolving with 
Glauber dynamics starting from a fully magnetized state, c): 
Qoo as a function of 9: the lines correspond respectively to 
Qoo(#) (solid) and Q^(9) (dotted) and the points correspond 
to the numerical values obtained for different nonequilibrium 
systems. The two values for 'Ising Id' correspond to the local 
magnetization (9 = 3/8) and to the global one (9 = 1/4). d): 
Same plot as in c) on a larger scale. 



ied the diffusing field </j(x, t) evolving according to the 
heat equation <9t<^(x, t) = V 2 y(x, t) with pbc in dimen- 
sion d starting from random initial condition (y(x, t — 
0)ip(x' , t = 0)) = <5 d (x — x'). It is known that the persis- 
tence exponent 9 = 9(d) associated to the diffusing field 
at the origin X(t) = <^(x = 0, t) depends continuously on 
d [H and in particular 9(d = 46) ~ 1 0. The proba- 
bility Q(t) can be easily computed numerically in any di- 
mension d by solving the heat equation and noticing that 
the field at the origin ip(0, t) can be simply written, for a 
large system size, as <p(0, t) ~ dr A d ~ 1 ^ 2 e~ r /*\&(r) 
where 'F(r) is a random field with short range correla- 
tions. We have computed Q(t) for d = 2,10,20,30,40 
and 50 and found, in all cases, Q(t) — > Qoo- The asymp- 
totic values Qoo as a function of 9, reported in Fig. [2t, 
are in good agreement, even quantitatively, with Q M (9) 
for multiplicative processes. In Fig. we have also re- 
ported the value of Qoo for another smooth process called 
the random acceleration process, for which 9 = 1/4. 
These data suggest that smooth processes, at variance 
with coarsening dynamics in Ising systems, are better 
approximated by multiplicative processes. 

A close look at Fig. [T] suggests investigation of other 
closely related cousins of l„ laK (t) defined in Eq. |T]), 
such as /x m ax(i) = max("ri,T2, • ■ ■ ,tn+i) or for instance 
Amax(i) = max(ri , T2 , • • • ,tjv). For renewal processes 
with Levy index 9 < 1, the analysis presented above 



can be extended to obtain exact results for the average 
of both observables. In both cases, the average grows 
linearly with t but with different ^-dependent prefactors. 
While the prefactor for the former case was computed in 
Ref. [l^] by a rather complicated but rigorous method, 
the latter case (A max (i)) has not been studied, to our 
knowledge, even for Brownian motion (9 = 1/2). We find 
(A m ax(i)} — Aoo((9)t which yields, in particular, a new 
constant Aoo(l/2) = 0.241749... for Brownian motion. 
The detailed studies of /i max (i) and A max (t) for other 
nonequilibrium processes will be reported elsewhere 0. 

In conclusion, we have shown that the average length 
of the longest excursion has rather rich and universal 
asymptotic time dependence for a variety of nonequi- 
librium processes. Our analytical and numerical results 
highlight the importance of extreme value statistics in 
generic nonequilibrium dynamics. 
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